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I review two subjects in which lattice simulations are making, or can make in the future, a significant 
contribution to particle physics phenomenology. The first subject is the evaluation of quantities which 
enter into the determination of the vertex A of the unitarity triangle from experimental measurements 
of decay rates and mixing amplitudes. These quantities include the form-factors for semileptonic B- 
decays, the leptonic decay constants of the B and Bg mesons and the S-parameters for B^-B^i and 
Bg-Bs mixing and the corresponding parameter for K-K mixing, Bk. In the second part of this talk I 
will review the status and prospects for the evaluation of if — » tttt amplitudes and for the subsequent 
study of the A/ = 1/2 rule and the evaluation of e'/e. 



1 Introduction 

In this lecture I will briefly review some of 
the contributions which lattice computations 
are making to particle physics phenomenol- 
ogy. The lattice formulation of quantum field 
theory together with large-scale numerical 
simulations is contributing to a wide variety 
of fundamental questions in particle physics, 
both theoretical and phenomenological. Here 
I will concentrate on one of the most impor- 
tant roles of lattice QCD, the evaluation of 
non-perturbative QCD effects in physical am- 
plitudes and other quantities. Indeed, it is 
frequently our inability to quantify the long- 
distance QCD effects in weak processes which 
is the dominant source of uncertainty in de- 
termining fundamental quantities from ex- 
perimental measurements and lattice simula- 
tions provide an ah initio framework for the 
evaluation of these effects. 

For many physical quantities lattice cal- 
culations have been performed for over ten 
years and the emphasis is on the reduction of 
systematic uncertainties. I will briefly outline 
the sources of some of the other uncertain- 
ties when presenting the results below, but 
let me now mention two important sources of 
error which the community is striving to re- 
duce. The first of these is quenching, the ne- 
glect of vacuum polarization and other quark 



loop-effects. Most large-scale phenomenolog- 
ical calculations have been performed in the 
quenched approximation, although increas- 
ingly calculations are being performed with 2 
flavours of sea quarks. Although it is natural 
to be skeptical about quenched calculations 
and consider full QCD (unquenched) ones as 
totally reliable, in my opinion neither of these 
views is fully justified at present. Where re- 
sults from quenched calculations can be com- 
pared with experimental measurements, they 
typically agree to within 10% or so, which 
is sufficiently accurate for some quantities 
and not so for others. On the other hand, 
it should also be remembered that in un- 
quenched simulations the masses of the sea 
and valence quarks are large (mjr/mp is typ- 
ically about 0.6 or more) so that significant 
extrapolations are needed, and this also leads 
to uncertainties. This brings me to the sec- 
ond source of lattice systematic error, which 
is increasingly being studied in detail, the ex- 
trapolation to the chiral limit. In order to 
avoid unphysical effects due to the finite vol- 
ume of the lattice, simulations are performed 
with up and down quarks with masses (m„ 
and md) in the region of rUg, the mass of the 
strange quark, and the results are then ex- 
trapolated to the physical values of m„ and 
rud (the computing cost also increases dra- 
matically as the masses of the quarks de- 
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crease) . For small values of the quark masses 
we can hope to exploit chiral symmetry (and 
chiral perturbation theory, in particular) to 
guide us in this extrapolation. The main 
question is therefore whether there is a region 
of overlap between the range of masses being 
used in simulations and those which are suf- 
ficiently light for chiral perturbation theory 
to apply. An added subtlety is that the chi- 
ral structure of the quenched theory is very 
different from full QCD. I will not describe 
these studies further, but they represent an 
important step towards the improvement of 
the reliability and precision of lattice compu- 
tations. 

I would like to stress that the material 
presented here represents only a small frac- 
tion of lattice results in general and lattice 
phenomenology in particular. The proceed- 
ings of the annual lattice conferences EJ con- 
tain detailed reviews of lattice contributions 
to different areas of phenomenology, as well 
as original contributions from the groups car- 
rying out the studies. For a summary of the 
future prospects in the subject I refer you 
to the report of an ECFA panel which was 
charged with the task of considering these □. 

In this talk I will focus on two topics, 
the contribution that lattice simulations are 
making to the determination of the vertex A 
of the unitarity triangle (section ^ and the 
status of lattice computations of the ampli- 
tudes for K — > TTTT decays and the prospect 
for the improvement in the precision of these 
calculations (section ||). 

The annual International Symposia on 
Lattice Field Theory provide an important 
forum for the collaborations to present their 
new results. This lecture was delivered before 
the 2001 Lattice conference (Latt2001), and 
in consequence the results presented in this 
written version were also largely compiled be- 
fore Latt2001. However, the new results pre- 
sented for kaon decays at Latt2001, and for 
e'/e in particular, are of considerable interest 
and I felt that it was necessary to present and 



discuss them. This is done in section 3.2 



2 Lattice QCD and the Unitarity 
Triangle 

In a number of talks at this conference we 
have seen the current status of the determi- 
nation of the vertex A of the unitarity trian- 
gle. Over-determination of the position of A 
is a convenient way of testing the consistency 
of the standard model of particle physics and 
of constraining its parameters. Theoretical 
inputs, and in particular quantitative esti- 
mates of non-perturbative QCD effects, are 
required in order to determine the possible lo- 
cus of the vertex from measurements of quan- 
tities such as the amplitudes of K'^ — K° mix- 
ing or studies of — mixing 0. Lattice 
QCD provides the opportunity for evaluat- 
ing these non-perturbative effects, and in ta- 
ble (taken from ref. i) I present a num- 
ber of the most important examples. The 
factors in bold-type in the second column of 
table 1^ are quantities which are frequently 
taken from lattice calculations. In the first 
half of this lecture I will review the status of 
the determination of these quantities. In ad- 
dition however, it should be noted that lattice 
calculations are also used in the determina- 
tion of the CKM matrix elements Vub and Vcb 
from exclusive semileptonic -B-decays, and I 
start with a brief review of these calculations. 

2.1 Exclusive Semi-Leptonic B -Decays 

The CKM matrix elements Vcb and Vub are 
determined from measurements of inclusive 
or exclusive semileptonic i3-decays. There 
have been a number of lattice computations 
of the exclusive form- factors, which combined 
with the experimental measurements of the 
amplitudes, allows the CKM matrix elements 
to be determined. The quark flow diagram 



"The determination of sin(2/3) from the mixing- 
induced CP-asymmetry in B — > J^/Ks is a beau- 
tiful and rare exception where there are essentially 
no hadronic uncertainties. 
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Table 1. Schematic table of some of the experimentpJ quantities which are measured from which information 
about the unitarity triangle is deduced (from ref. 13). Factors in bold- type represent quantities which are 
calculated in lattice computations. 
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Figure 1. Schematic representation of semileptonic 
S-decays. 
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for exclusive semi-leptonic B decays is rep- 
resented in fig. |l|. Lorentz and parity invari- 
ance allows us to write the decay amplitudes 
in terms of invariant form factors, e.g. for a 
decay into a pseudoscalar P {P — D or ir for 
example) we can write: 



— 

{P{pp)\V,mB{pB)) = fo(q^) ^ 



+f+(q') 



Ml - Ml 



9m 



, (1) 



where in this case there are two form-factors 
/o and /_!-. Parity invariance implies that 
only the vector component V of the weak 
V ~ A current contributes when the final- 
state hadron is a pseudoscalar. For B vec- 
tor decays, both the vector and axial-vector 
currents contribute and there are four form- 
factors {Ai,2,A{= Ao-A3),V). 

B D''*^ Iv- Decays: Lattice calcu- 
lations could, in principle and perhaps also 
in practice, make a contribution to the de- 
termination of Vcb by determining the corre- 
sponding form factors. This is a challenging 
task however, since, in order to make an im- 
pact, one needs to calculate small corrections 
to the result in the heavy quark limit. In 



particular, for B p decays we have: 
dT{B D*lv) 



= ^x\VcbVT\u:) , (2) 



where K is a known kinematic factor, uj = 
vb ■ vjj- {vb and are the four- velocities 
of the B and D* mesons respectively) and 
= 1-1- corrections. In order to make a 
contribution one has to be able to calculate 
the 1 / itlq non-perturbative corrections to the 
distribution at zerojecoil (lj = 1). There has 
been a suggestion □ that by calculating the 
ratio of ratios 

{D\c^^b\B) {B\hoc\D) 
{D\cjoc\D) {B\hob\B) ^ ' 

it may be possible to determine the form- 
factors sufficiently accurately. This clearly 
requires an excellent control of the system- 
atic errors and it remains to be seen whether 
sufficient precision will be possible. I men- 
tion in passing that lattice simulations have 
been performed which reproduce the CLEO 
data on the oj distribution I will however, 
focus more on the b ^ u decays where the 
opportunity for lattice calculations to make 
a contribution to phcnomenological studies is 
greater. 

B piv-Decays: Simulations of i3 — > tt 
or p semileptonic decays yield the form fac- 
tors at large values of . The reason for 
this limitation is that the momenta of the 
pion or p-meson must be kept small in or- 
der to avoid artefacts due to the granular- 
ity of the lattice. The calculations have been 
performed for many years now, and in fig. |^ 
I reproduce a figure from Claude Bernard's 
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Figure 2. Recent data for the form- factors /+ and fo 
of semileptonic B — > irlu decays from C. Bernard's 
review talk at Lattice 2000 □. 
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review talk at last year's lattice conference, 
showing results for the B tt form fac- 
tors from four collaborations in the region 
19GeV2 <q^ < 23GeV2. The collaborations 
use different formulations for the &-quark, but 
the results are in reasonable agreement. 

Although much effort is being devoted 
to extrapolating the lattice results to smaller 
values of q^, using as many theoretical con- 
straints as possible (such as heavy-quark 
symmetries, unitarity and analyticity, kine- 
matical constraints and soft-pion relations), 
the most meaningful applications of lattice 
results are (and/or will be) to the experimen- 
tal distributions at large values of g^. A pre- 
liminary comparison of the UKQCD lattice 
data to the CLEO measurement of the con- 
tribution to the width for B p decays from 
the region lAGeV'^ < < 20.3 GeV^ (which 
corresponds to CLEO's large q^ bin) gives: 



Ar = (7.ill:g)|K6rio 

UKQCD \ (4) 
= (7.1 ±2.4) X 10^ s"^ 

CLEO i, (5) 

from which one obtains |Kif,| — (3.2 ± 0.6) x 
10"'^. As the experimental statistics in- 
creases and lattice results get more precise, 
such comparisons will become the most di- 
rect way of determining Vub from exclusive 



Figure 3. Schematic representation of the leptonic 
decay of the B-meson. 
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B 



decays. 

2.2 The Decay Constants Jb and fs^ 

In fig. ^ I show the quark-flow diagram for 
the leptonic decay of the _B-meson. The non- 
perturbative QCD effects are contained in the 
matrix element 



[0\A^{0)\B{p))^^fBP^.^ 



(6) 



where is the axial-vector current with 
the appropriate flavour quantum numbers. 
Using Lorentz and Parity Invariance we see 
that all the nonperturbativc QCD effects are 
parametrized in terms of a single number 
Jb the (leptonic) decay constant of the B- 
meson ^. Quenched calculations of /b have 
been performed for about 15 years now and 
a careful analysis of all the systematic errors 
(apart from quenching) is possible (see for ex- 
ample C. Bernard's review at Lattice 2000 □). 
In fig. ^ I update Claude Bernard's compila- 
tion of recent results using a variety of differ- 
ent formulations for the b quark His con- 
clusion from these results for the quenched 
value of fs (with which I concur) is: 

/s,quenchcd = 175±20MeV, (7) 

and this is shown as the shaded box in fig. |3[ 
It may appear a little strange that with 
so many results being included in fig. ^, the 
final error is as large as in eq. (Q) and as indi- 
cated in the shaded region of the figure. This 
is a manifestation of the difficulty in con- 
trolling the systematic uncertainties, many of 
which are common to the different determi- 
nations and the estimate of the overall error 



''The convention for the normahzation used here cor- 
responds to /tt — 132 MeV. 
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Figure 4. Results for fg from various groups in the 
quenched approximation a. Statistical and system- 
atic errors have been combined in quadrature. 
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requires a careful analysis of the treatment of 
these uncertainties by each group . For ex- 
ample, we know that quenching induces er- 
rors of 0(10%) in many physical quantities. 
In particular the value of the lattice spacing 
determined by using different physical quan- 
tities to set the scale (e.g. mp or /jr) typi- 
cally also varies by this amount. Thus there 
is an irreducible error in the value of fs in 
the quenched approximation of about 10% 
(or about 20 MeV). 

The emphasis is now turning to un- 
quenched calculations, see fig.|[ There is 
some belief that fB,Nj=2 is 10-15% larger 
that the decay constant in the quenched ap- 
proximation and C. Bernard's conclusion is 



Figure 5. Results for /n-from 4 simulations using two 
flavours of sea quarks l3. Statistical and systematic 
errors have been combined in quadrature. 
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An important parameter which appears 
in analyses of the Unitarity Triangle is 



(9) 



where the B-parameters of neutral i?-meson 
mixing are defined in section 2.3 below. It 



will be shown below that lattice simulations 
indicate that the B-parameter varies slowly 
with the light-quark mass and so it is inter- 
esting and instructive to consider the ratio 
/Bs / fsd ■ The lattice results for this quantity 
have been very stable and C. Bernard con- 
cludes from these results that 

JB J quenched 
Ib 

in the quenched theory and in full QCD re- 
spectively. 

Thus the ratio in eqs. ( |l0|) and (|ll| ) is de- 
termined rather precisely. This should not be 
too much of a surprise since the key point to 
note is that for both ^ and fsJlBa it is the 
difference from 1 which is being computed. 



1.15 ±0.04 (10) 

1.16 ±0.04. (11) 
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Figure 6. Schematic representation of B^-B" mixing. 
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Lattice errors of 30% (which is a conserva- 
tive estimate) therefore correspond to errors 
of only about 5% on ^. 

In this short lecture I do not have the 
opportunity of explaining the formulations of 
heavy quarks on the lattice in any detail, but 
let me mention briefly the reason why differ- 
ent groups use different formulations. The 
number of lattice points is limited by the 
available computing resources and we require 
the volume of the lattice to be larger than the 
hadrons being studied. We also require the 
lattice spacing a to be small enough to avoid 
discretisation errors (i.e. artefacts due to the 
granularity of the lattice) so that the choice of 
the value of a is a compromise between two 
sources of possible error. Typically, in cur- 
rent simulations, one takes ^ 2-3GeV. 
This is less than the mass of the 6-quark, nib, 
so that we cannot study the propagation of a 
physical &-quark on presently available lat- 
tices. To circumvent this difficulty, lattice 
results for 6-physics are obtained by taking 
those obtained with the heavy-quark mass, 
mq, in the region of the mass of the charm 
quark (toq ~ rric) and performing the extrap- 
olation niQ —^ nih or by performing simula- 
tions in effective theories, such as the Heavy 
Quark Effective Theory or Non-Relativistic 
QCD or by a combination of these two ap- 
proaches. There is therefore reasonable con- 
trol over this source of systematic uncertainty 
and as computing resources increase we will 
be able to verify explicitly that this is indeed 
the case. 



2.3 B-B Mixing 

The quark flow diagram for B-B mixing is 
drawn in fig.^. The non-perturbative QCD 
effects in this important process are con- 
tained in the single matrix element: 

M{p) = ( 5° I 67^(1-75)9 67^(1-75)9 \Bl), 

(12) 

where q represents the d or s quark. It 
is convenient and conventional to define B- 
parameters by factorizing the vacuum satu- 
ration contribution: 

M{^,)^plm%BBM- (13) 

The Bb^ (m) parameters are renormalization 
scheme and scale-dependent and therefore it 
is again convenient and conventional to define 
scheme-independent (up to NLO) quantities 



rynlo 



"s(/-i) 



2//3o 



1 + ^^J„, SB,(/i), 

(14) 

where Jn^ is a known constant calculated in 
perturbation theory. 

Little has changed since last year, when 
C. Bernard's summary of the corresijonding 
results at the lattice conference was □: 



Bb^ = 1.30 ±0.12 ±0.13 



fB,\/BB, = 230±40MeV 



Bf 



1.00 ±0.04 



^ = 1.16 ±0.05 



(15) 
(16) 

(17) 

(18) 



2.4 K°-K° Mixing and Bk 



The quark flow diagram is similar to that for 
Bb above (see fig. |^). The non-perturbative 
QCD effects are contained in the matrix ele- 
ment: 

(i^°|(s-7;.(l-75)rf)(s7p(l-75)d)|i^°) . (19) 

Chiral symmetry plays a central role in 
the determination of the corresponding B- 
parameter Bk- This presents a relative dif- 
ficulty for simulations performed using the 
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Wilson formulation of lattice fermions (or ex- 
tensions of this formulation), since in this 
case we don't have explicit chiral symmetry 
until the extrapolation to zero quark mass is 
performed. Specifically, the operator in the 
matrix element of eq. (jl^) mixes with other 
operators of dimension 6 and the matrix el- 
ements of these operators have to be sub- 
tracted, leading to a loss of precision. Over 
the years however, techniques have been de- 
veloped to perform these subtractions non- 
perturbatively and reasonably effectively. 

There are two recent and related pro- 
posals to circumvent the need for the sub- 
traction of the additional operators, based 
on Bernard's observation that CPS-symm- 
etry (where S=s <-!■ d, the interchange of the 
s and d quarks) implies that the parity-odd 
component of the AS — 2 operator in eq. ( p^ 
renormalizes multiplicativcly. However, it is 
the parity-even component which is non-zero 
in eq. (|l|). 

The first proposal is to use twisted-mass 
QCD 0: 

(20) 

where "0 represents the isodoublet of light- 
quarks (r^ is a matrix in this space), 
is a parameter and Dw is the Wilson for- 
mulation of the Dirac operator. The (mul- 
tiplicatively renormalized) parity-odd opera- 
tor in this theory corresponds to the physical 
(parity-even) operator in QCD. 

It is also possible to use a (chiral) Ward 
Identity to determine the physical matrix el- 
ement of the parity-even AS = 2 operator 
from a measurement of that of the parity- 
odd component, Oy^^^, without the twisted 
mass-term El. A recent simulation with this 
formulation gives: 



BK(2GeV) = 0.73 ±0.07 



+0.05 
0.01- 



(21) 



Nevertheless, the most precise evaluation 
of Bx comes from simulations using the stag- 
gered formulation of QCD, in which chiral 
symmetry is explicit at the expense of a more 



Figure 7. Results for B^^^{2GeY) in the con- 
tinuum limit from variopi groups (mainly) in the 
quenched approximation Oj. Statistical and system- 
atic errors have been combined in quadrature (where 
these have been presented). 
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complicated flavour structure. The evalua- 
tion of Bk is also an excellent testing ground 
for new formulations of chiral fermions based 
on the Ginsparg- Wilson relation and results 
of Bk = 0.787(8) and 0.737(11) have_been 
reported by the CP-PACS M and RBC Q col- 
laborations using Domain Wall Fermions. 

Fig. contains a compendium of recent 
results for Bk compiled by L.Lellouch, the 
rapporteur at 2000 Lattice conference, whose 
conclusion from these results is c3: 

^NDR(^2 GeV) = 0.628 ± 0.042 ± 0.099 



B 



NLO 

K 



= 0.86 ±0.06 ±0.14 . (22) 



The errors include estimates of those due to 
quenching obtained using Sharpe's analysis 
with quenched chiral perturbation theory EH. 

In extensions of the standard model, 
and in particular in supersymmetric theo- 
ries, AS* — 2 operators other than that in 
eq. ( |l9[ ) also contribute to Bk- These can 
also be evaluated in lattice simulations Ei'0 
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thus helping to provide constraints on the 
properties of supersymmetric models. 

3 K ^ TTTT Decays 

The evaluation of the amplitudes for nonlep- 
tonic weak decays, particularly of B and K 
mesons, represents a major challenge for lat- 
tice physicists. At this conference we have 
heard many interesting experimental results 
for two-body exclusive decays of B-mes ons 
such as i? — > TTTT or i? ^ nK decays c3, 
but, at present, we are unable to perform 
lattice calculations of the corresponding ma- 
trix elements. Considerable progress is being 
made however, towards reliable calculations 
of K TTTT decays, including quantitative 
studies of the A/ — 1/2 rule and an evalu- 
ation of e'/e. In this section I will outline 
some of this progress and attempt to convey 
the optimism we feel for future prospects. For 
a more comprehensive and detailed recent re- 
view see ref. Ell and references therein. 

The lattice contribution to the evalua- 
tion of K ^ TTTT amplitudes begins with the 
use of the operator product expansion lead- 
ing to an expression for the AS' = 1 effective 
weak Hamiltonian in terms of Wilson coefR- 
cient functions Ci and renormalized local op- 
erators Oi{fi): 



n,AS= 



G 



i = --|5]a(A.)0.(/.). (23) 



The CiS contain the perturbative QCD ef- 
fects which give the evolution from the mass 
of the W to the (perturbative) renormaliza- 
tion scale /x. The non-perturbative physics is 
contained in the matrix elements (7r7r|Oi|isr), 
and the role of lattice simulations is to eval- 
uate these matrix elements. 

From lattice computations one obtains 
the matrix elements of the bare lattice oper- 
ators with the lattice spacing a as the ultra- 
violet cut-off. From these we must construct 
the finite matrix elements of renormalized op- 
erators, and this ultra-violet problem is, in 



principle at least, fully solved E3. Several 
non-perturbative techniques have been devel- 
oped to determine the cqpre^ponding renor- 
malization coefficients ocJEj (for recently 
calculated renormalization constants in per- 
turbation theory for the Domain Wall formu- 
lation of lattice fermions see ref. E£l). 

Two main approaches are used to deter- 
mine the decay amplitudes from lattice sim- 
ulations: 

i) The K ^ n (and K — > vacuum) ma- 
trix elements are computed directly, and the 
K — > TTTT matrix elements are obtained using 
soft-pion theorems and (lowest order) chiral 
perturbation theory. 

ii) The K — > tttt matrix elements are com- 
puted directly. 

Although, at first sight, it may appear that 
the second approach is clearly better (and 
I expect that it will eventually become the 
standard one), it does involve a two-hadrou 
final state which presents some subtleties cZl. 
We are trying to determine physical decay 
amplitudes from matrix elements computed 
in a finite Euclidean volume. Lellouch and 
Liischer have initiated substantial progress 
towards the solution of this infrared problem 



and I will briefly review this in section 3.1 . In 
section 3.2 I will review some recent numer- 



ical results and discuss prospects for future 
calculations of K ^ tttt amplitudes. 

3.1 K TTTT Decays in a Finite Volume 

The infrared problem for the evaluation of 
K — > TTTT decay amplitudes arises from two 
sources, the unavoidable continuation of the 
theory to Euclidean space-time (the Maiani- 
Testa Theorem E3) and the use of a finite 
volume in numerical simulations. An impor- 
tant step towards the solution of this problem 
has been achieved by Lellouch and Liischer @ 
who derived a relation between the K — > tttt 
matrix element in a finite volume and the 
physical decay amplitude: 

\{nn\nw{o)\K)\'' = \v{nn\nw{Q)\K)v? 
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In the first line of eq. (|2 



{qcf>'{q) + k5'o{k)}k=k^.{2A) 



the left-hand 
side is the infinite-volume matrix elements 
and the right-hand side is the finite-volume 
(which might, for example, be computed in 
a lattice simulation). The second line is the 
factor relating these and the main message 
that I am trying to convey here is that there 
is a known factor which relates the physical 
amplitude and the finite-volume matrix ele- 



ment. In eq. (|24| ) k is related to the centre of 
mass energy W by: 

W = 2-y/mf~4~A?2 and fc^r = ^^Jrr^^~im^ ; 

(25) 

q — kL/2'K where L is the length of the lat- 
tice, (/)(g) is a known kinematic function of 
q and is a consequence of the cubic shape 
of the finite lattice and 5{k) is the physical 
(infinite-volume) s-wave phase-shift (the ex- 
plicit formula in eq. (^4|) applies only to this 
partial wave). The remaining finite- volume 
corrections decrease exponentially as the vol- 
ume increases. 

In lattice simulations we calculate cor- 
relation functions at large (Euclidean) times 
so as to isolate the ground state. Most fre- 
quently, as for example in the computations 
used to obtain the results presented in sec- 
tion 1^, the interpolating operators used in 
the correlation functions are such that the 
ground state corresponds to the lightest par- 
ticle with some specified quantum numbers. 
For K — > TTTT decays we are also interested in 
two-particle states. One of th£^ consequences 
of the Maini- Testa theorem EZl is that such 
correlation functions are dominated at large 
times by the states in which the kaon and 
each of the two-pions arc (almost) at rest. 
Such a kinematical situation is clearly un- 
physical. Lellouch and Liischer §1 make the 
interesting and significant observation that, 
since energy levels in a finite box are dis- 
crete, it is possible to tune the size of the 
box (lattice) in such a way that the energy 



of the first excited two-pion state is precisely 
the mass of the kaon. This would require a 
volume of about 6 fm, somewhat larger than 
currently used, but one which should become 
accessible in quenched simulations with the 
next generation of dedicated computers. Al- 
though the first excited state is more difficult 
to extract than the ground state, this is not 
likely to present a major difficulty. 

Lellouch and Liischer derive the formula 
in eq. (|4|) for a lattice with a fixed large 
volume V , chosen in such a way that the 
decay of the kaon with physical kinematics 
corresponds to one of the two-pion energy 
levels accessible on this lattice (such as the 
first excited state mentioned above). We 
have recently rederived eq. ([24[), taking the 
— > cxD limit at fixed physics starting from 
the Liischer quantization condition ?1 



(t>{<l) +5o{k) = nvr. 



(26) 



which gives the spectrum of _Lwo-particle 
states in a finite cubic volume t2l. We were 
able to establish the validity of eq. ( ]2^ ) for 
all elastic states below the inelastic thresh- 
old, with exponential accuracy in_the volume, 
extending the derivation in ref.E3 which was 
presented for the lowest seven levels. We were 
also able to demonstrate that the formula 
was valid when the two-pion energy does not 
match the mass of the kaon ttlk and the in- 
serted operator (e.g. one of the Oi in the op- 
erator product expansion for the weak Hamil- 
tonian eq. ( ^3|) ) carries non-zero energy and 
momentum. This is particularly useful in at- 
tempts to determine the coefficients of the 
operators appearing in higher orders of the 
chiral expansion (see sec. |3.2| below). 

3.2 Status and Prospects for the 
Evaluation of K tttt Decays. 

In this section I present some recent lat- 
tice results of the matrix elements which 
contribute to i^T nn decay amplitudes 
and discuss the exciting prospects for further 
progress in improving the precision. Since 
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Figure 8. Two contributions to nonleptonic kaon de- 
cay amplitudes. 




a) Disconnected Emission 




b) Disconnected Penguin 

the delivery of this lecture in July, two groups 
have presented results which include the eval- 
uation of A/ = 1/2 amplitudes including the 
determination of e'/e, and in view of the in- 
terest in such calculations I will also comment 
briefly on these results. 

A number of different mechanisms con- 
tribute to the amplitudes for K ^ tt-k decays, 
and in fig. H I show some examples. The four- 
quark vertex at which there is an s — s- u or 
s ^ d transition represents the insertion of 
one of the four-quark operators appearing in 
the effective Hamiltonian in eq. (^) . Lattice 
calculations have shown that it is not possi- 
ble to explain the A/ — 1/2 rule, i.e. the 
experimentally observed enhancement of the 
amplitudes for A/=l/2 decays by a factor of 
about 22 relative to those for A/=3/2 de- 
cays, with emission diagrams only. I start 



with a discussion of attempts to evaluate the 
penguin contractions using Wilson-like lat- 
tice fermions. In order to obtain the phys- 
ical contribution from the penguin diagrams, 
in general we have to subtract large unphys- 
ical artefacts, terms which diverge as inverse 
powers of the lattice spacing (power diver- 
gences). From the diagram in fig. ||(b) we 
can see that the inserted four-quark operator 
can mix with quark bilinears of the type dTs 
(where T is one of the Dirac matrices) unless 
there are symmetries to prevent this mixing. 
Although lattice symmetries do soften the di- 
vergences corresponding to the mixing, large 
subtractions (direct or indirect) are neverthe- 
less unavoidable. This is the reason for the 
absence up to now of sufficiently precise re- 
sults for A/ — 1/2 decays. 

The excitement in the SPQqqR (South- 
ampton-Rome- (QCD)-Paris) collaboration, 
of which I am a member, is due to the fact 
that for the first time we have a lattice sig- 
nal for the amplitudes, which we will analyse 
to determine e' /e and study the A/ = 1/2 
rule. This is partly due to improved theoreti- 
cal techniques to reduce the subtractions and 
to deal with them non-perturbatively, and 
partly due to improved computing facilities. 
In fig. H I show the (preliminary) raw lattice 
data for the ratio of correlation functions as 
a function of the time t, from which one of 
the relevant matrix elements, 7r7r(7r7r|0_ jif), 
is determined. There is a stable region in t, 
where the matrix element can be seen to be 
non-zero and we will increase the statistics to 
reduce the error which is currently still large. 



In fig. 10 I illustrate the huge subtrac- 



tions which generally have to be performed 
with data from the CP-PACS collaboration 
for the matrix element of Og, which is one of 
the important operators contributing to e'/e. 

I now turn to the matrix elements of 
A/ — 3/2 operators between the kaon and 
two-pion states. These can can now be eval- 
uated with good precision, and as an example 
I present in fig. |l^ the A/ — 3/2 matrix ele- 
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Figure 9. Raw lattice data for the matrix element 
i—o{ttit\0~ \K) as a function of the time. 



ments of the electroweak operator Og which 
is one of the key components in the evalua- 
tion of e'/e. We are currently undertaking a 
detailed study, up to next-to-leading order in 
the chiral expansion for the AI = 3/2 matrix 
elements of these operators (04.7,8)- Prelimi- 
nary results indicate that it is possible to ob- 
tain a AI = 3/2 if —^ TTTT decay amplitude, 
which is nevertheless consistent with a large 



Bji (see sec.2.4). 

There is a long history of lattice stud- 
ies of kaon decays in which operator matrix 
elements of the type {tt\0\K) are computed 
and combined with soft-pion theorems and 
chiral perturbation theory to obtain the de- 
cay amplitudes. Here I mention one such re- 
cent result, for the matrix element of the elec- 
troweak penguin operator Og 113, for which 
lattice results give significantly smaller val- 
ues than other determinations, e.g. 

\i=2{TnT\Os\K°)\ = (0.5±0.1)GeV3 , (27) 

in the NDR renormalization scheme at a scale 
of 2 GeV (such a small result was also con- 
firmed from computations of K — > tttt ampli- 
tudes by the SPQ^^j^R collaboration in the 
presentation of their preliminary data at Lat- 
tice 2001 @). This can be compared to the 
larger values, e^. 2.22±0.67GeV3 El and 
(3.5±1.1) GeV^ Efl obtained using other (con- 



Figure 10. An example of the very large subtrac- 
tions present in some lattice evaluations of the AI = 
1/2 transitions. The figure shows the values of the 
unsubtracted and subtracted bare matrix elements 
j_o(7r7r|06 \K) as a function of the quark massifrom 
a simulation by the CP-PACS Collaboration EiJ. 
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tinuum) techniques such as dispersion rela- 
tions and the 1 /Nc expansion ^. There is also 
a recent continuum result of 1.2±0.5Gev^ 
obtained using spectral functions. A large 
value of the matrix element of Og would re- 
quire very large values indeed for the matrix 
elements of Og in order to explain the mea- 
sured value of e'/e and so it will be interesting 
to follow the developments in these calcula- 
tions. 

As mentioned above, two groups have re- 
cently presented results for e'/e and other 
properties of K tttt decay amplitudes, 
from a computation of K t: matrix ele- 
ments using the Domain Wall Fermion for- 
mulation for lattice quarks (these results were 
presented since this lecture was delivered). 
In view of the experimental result E e'/e = 
(17.2±1.8) lO-'^ these results may seem some- 
what disturbing: 



e'/e = (-8-^-4) 10"^ 
e'/e = (-7^-2)10-^ 



RBC38 (28) 
CP-PACsi].(29) 



After irteliuding the recently cateulated O(a^) cor- 
rections the authors (ji^ef. Eil find their results 
becomes (2.4 ± 0.8) Gev^ M. 
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Figure 11. Preliminary data for the matrix element of 
the electroweak penguin operator Og for a particular 
choice lei the quark mass as a function of the lattice 
time t E3. 
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Of course, it would be very exciting to be 
able to confidently deduce the existence of 
new physics from the discrepancy between 
the lattice results in eqs. (|2|) and (|2|) and 
the experimentally measured value of e'/e. 
Such a conclusion remains a tantalizing pos- 
sibility. Before this can be done however, we 
need to be reassured that the systeniatics of 
the computations are sufficiently under coiLj 
trol.^lthough the results for e'/e from refs.Ea 
and El are consistent with each other, this is 
not the case for some other quantities, includ- 
ing some of the separate components in e'/e. 
For example both groups find a large, but 
different, value for the ratio Re(ylo)/Re(j42), 
where and 2 denote the isospin of the two- 
pion system (and hence an enhancement of 
A/ = 1/2 decays!). The RBC collabora- 
tion find ReAo/ReA2 ^24-27 (remarkably 
close to the experimental value) whilst CP- 
PACS find ReAo/ReA2 =9-^-12. Since both 
groups use similar methods (but with some 
important differences, for example in the nor- 
malisation of the operators) we need to un- 
derstand the reason for discrepancies such as 
these. 

It has to be stressed that this calcula- 
tion is much more difficult and subtle than 



those reported in section Q. Given the huge 
subtraction of power divergences, illustrated 
in fig. |l^, a good understanding of the chi- 
ral properties of the theory with the param- 
eters used in the simulation is crucial. The 
computations were performed using Domain 
Wall Fermions with N^^ — 16 points (where 
iVs is the number of points in the fifth direc- 
tion), and although both groups claim that 
this is sufficient for the residual chiral symme- 
try breaking effects to be fully under control, 
it would be very reassuring to confirm this 
by increasing whilst keeping the other pa- 
rameters fixed. We also need to be able to un- 
derstand the consequences for these calcula- 
tions of the recent observation by Golterman 
and Pallante cB that in the quenched approx- 
imation there are additional (spurious) chi- 
ral logarithms. It should also be remembered 
that the results for K — )■ tttt amplitudes in 
eqs. eqs. (^8|) and ( ^ ) were obtained from 
the determination of — > tt matrix elements 
using lowest order chiral perturbation theory, 
and one can ask about the precision of this 
procedure. Nevertheless, in spite of caveats 
such as these, these new results are very ex- 
citing and mark the beginning of a new era 
in lattice studies of kaon decays. 

4 Conclusion 

Lattice simulations provide the opportunity 
to evaluate non-perturbative QCD effects 
from first principles with no model assump- 
tions or parameters. There is a large range 
of quantities of central importance to parti- 
cle physics which are being computed in lat- 
tice simulations (indeed it is far too large a 
range to be considered in a review talk of 30 
minutes). For some quantities, such as those 
which enter into the analysis of the Unitarity 
Triangle which were discussed in section |2[ 
the emphasis is now on the reduction of sys- 
tematic errors. For others, such as the eval- 
uation of nonleptonic weak decays in general 
and isT — > TTTT decays in particular, we are 
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still learning how best to extract the physical 
quantities. The range of quantities which can 
be studied is constantly expanding. 
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